Draft Version September 21, 2012 

Preprint typeset using WT^K. style emulateapj v. 5/2/11 



o 

CO 



o 

u 

CM 

o 

>: 

00, 
OS' 

in 

: 

00' 

o: 



SPATIAL ANIS0TR0PY OF GALAXY KINEMATICS IN SLOAN DIGITAL SKY SURVEY GALAXY 

CLUSTERS 

Andreas Skielboe 1 ' 2 , Radoslaw Wojtak 1 , Kristian Pedersen 1 ' 2 , Eduardo Rozo 3 ' 4 ' 7 , and Eli S. Rykoff 5 ' 6 

Dark Cosmology Centre, Niels Bohr Institute, University of Copenhagen, DK-2100 Copenhagen, Denmark 
Space Science Center, Niels Bohr Institute, University of Copenhagen, DK-2100 Copenhagen, Denmark 
Department of Astronomy and Astrophysics, The University of Chicago, Chicago, IL 60637, USA 
Kavli Institute for Cosmological Physics, The University of Chicago, Chicago, IL 60637, USA 
SLAC National Accelerator Laboratory, Menlo Park, CA 94025, USA and 
Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA 
Draft Version September 21, 2012 

ABSTRACT 

Measurements of galaxy cluster kinematics are important in understanding the dynamical state 
and evolution of clusters of galaxies, as well as constraining cosmological models. While it is well 
established that clusters exhibit non-spherical geometries, evident in the distribution of galaxies on 
the sky, azimuthal variations of galaxy kinematics within clusters have yet to be observed. Here we 
measure the azimuthal dependence of the line-of-sight velocity dispersion profile in a stacked sample 
of 1743 galaxy clusters from the Sloan Digital Sky Survey (SDSS). The clusters are drawn from the 
SDSS DR8 redMaPPer catalog. We find that the line-of-sight velocity dispersion of galaxies lying 
along the major axis of the central galaxy is larger than those that lie along the minor axis. This is 
the first observational detection of anisotropic kinematics of galaxies in clusters. We show that the 
result is consistent with predictions from numerical simulations. Furthermore we find that the degree 
of projected anisotropy is strongly dependent on the line-of-sight orientation of the galaxy cluster, 
opening new possibilities for assessing systematics in optical cluster finding. 
Subject headings: galaxies: clusters: general — galaxies: kinematics and dynamics 



1. INTRODUCTION 

Clusters of galaxies constitute the largest, gravitation- 
ally collapsed, structures in the universe. They offer a 
unique opportunity to study the formation and evolu- 
tion of structure on cosmological scales. Additionally, 
the mass distribution of galaxy clusters can be obser- 
vationally probed both dynamically and via weak lens- 
ing, making galaxy clusters ideal laboratories for study- 
ing d ark matter and modificat ions to general relativity 
(e.g. Rapetti et al.l (|2010l 12012ft ). 

Observations and simulations clearly show that 
clusters exhibit triaxi al rather than spherical 
shapes dBinggelil 119821 : iBailin fc SteinmetzT 120051 
iGottlober fc Yepesl 12007^ iHoekstra et al.l (|2004ft made 
the first detection of the flattening of galactic dark 
matter halos using weak lensing and ? measured the 
azimuthal variation of the weak gravitational lensing 
signal around galaxies, constraining galactic dark matter 
halo ellipticity. Non-zero ellipticity of dark matter halos 
has also been confirmed in strong lensing measurements 
by studying the angular distribution of giant arcs 
around groups and c lusters of galaxies ([Dalai et al.l 
l200l iMore etaIll20T2h . 

The velocity fields of galaxies in clusters are known 
to feature non- trivial radial dependencies, but azimuthal 
variations are often overlooked. If the shape of the under- 
lying dark matter halo is indeed multiaxial, it is natural 
to think that it will be reflected in an anisotropic velocity 
field. Indeed this has b een confirmed and studied using 
numerical simulations (jTormenl IT 9 9 71; iKasun fc Evrardl 
2005: IWhite et aLll2Q10t iSaro e t al. 2012ft . Because clus- 
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ter mass estimators often assume spherically symmetric 
velocity fields, characterizing the degree of anisotropy in 
galaxy clusters is crucial in understanding its impact on 
cluster mass estimates. It may also provide clues to the 
ongoing formation of galaxy clusters. 

The effect of an anisotropic velocity field should man- 
ifest itself in an azimuthal variation of the projected ve- 
locity dispersion of member galaxies. To our knowledge 
no attempts have been made to measure the azimuthal 
variation of galaxy velocity dispersions in a large sample 
of galaxy clusters. 

The aim of this Letter is to test for azimuthal de- 
pendence of the projected velocity dispersion of cluster 
galaxies in a stacked sample of galaxy clu sters. We use 
the S loan Digital Sky Survey (SDSS) DR8 (jAihara et al.l 
12011ft redMaPPer cluster catalog (described below) to 
search for an anisotropic velocity field by measuring the 
velocity dispersion of cluster galaxies along the major 
axis of the central galaxy (CG), and galaxies along the 
corresponding orthogonal axis. The results are compared 
with expectations from numerical simulations. We adopt 
a flat ACDM cosmology with Q rn = 0.3 and Hq — 70 km 
s" 1 Mpc -1 . 

2. METHODS 

The expected degree of the relative difference in dis- 
persion is of the order of a few percent (see Section I3TT]) . 
The strategy for detecting a signal of this magnitude is 
to stack a large number of galaxy clusters. 

2.1. The redMaPPer Cluster Catalog 

Our study utilizes the SDSS DR8 red-sequence 
matched- filter probabilistic percolation ("redMaPPer") 
cluster catalog (E. S. Rykoff et al. 2012, in preparation), 
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based on the optimized r ichness estimator A ([Rozo et al.l 
120101 : iRykoff et al1l2Q12h . redMaPPer is a photometric 
cluster catalog that identifies galaxy clusters as over- 
densities of red-sequence galaxies. The algorithm is di- 
vided into two stages: a calibration stage, where the red- 
sequence model is determined directly from the data, and 
a cluster-finding stage. These two stages are iterated sev- 
eral times before a final cluster-finding run is made. 

In the calibration phase, redMaPPer empirically cal- 
ibrates the color distribution (mean and scatter) of 
red-sequence galaxies as a function of redshift and 
magnitude. The calibration itself is also an itera- 
tive procedure. First, we use red galaxies from the 
SDSS lu minous red galaxy jEisenstei n et al.l [2001) and 
"main" (Str auss et al.l 120021 ) spectroscopic galaxy sam- 
ples (over 50% of the sky) as training galaxies. These 
galaxies are used as "seeds" to look for significant over- 
densities of nearby galaxies with the same color as the 
seed galaxy. The resulting cluster list is used to fit a full 
red-sequence model including zero point, tilt, and scat- 
ter. Moreover, the algorithm utilizes all colors (u — g, 
g — r, r — i, and i — z) simultaneously, so the "scatter" 
is characterized not by a single number but by a covari- 
ance matrix. The red sequence model is calibrated down 
to a luminosity threshold of 0.2 L* at the cluster red- 
shift (the optimal depth for cluster richness estimation 
IRykoff et al.l l2Q12h . thereby leveraging the bright spec- 
troscopic galaxy sample to obtain a model of the red 
sequence that extends to dim magnitudes. 

With our model of the red sequence at hand, we pro- 
ceed to run the cluster-finding algorithm as described 
below. The resulting cluster catalog is then used to re- 
calibrate the model of the red-sequence as above, and 
the procedure is iterated until convergence in the red se- 
quence model parameters is achieved. At that point, we 
perform one final cluster-finding run to arrive at the final 
redMaPPer catalog. 

Given a red-sequence model, the cluster finding pro- 
ceeds as follows. First, we consider all photometric 
galaxies as candidate cluster centers. We use our red- 
sequence model to calculate a photometric redshift for 
each galaxy, and evaluate the goodness of fit of our red 
galaxy template. Galaxies that are not a reasonable fit 
to the red-sequence model at any redshift are immedi- 
ately discarded. For the remaining galaxies, we use this 
initial redshift guess and the red-sequence model to eval- 
uate the richness (A) and cluster like lihood using a mult i- 
color generalization of the method in lRvkoff et al.l ((2012). 
When a significant number of red galaxies are detected, 
we re-estimate the cluster redshift by performing a simul- 
taneous fit of all the high-probability cluster members 
to the red-sequence model. This procedure is iterated 
until convergence is achieved between member selection 
and cluster photometric redshift. The resulting list of 
candidate cluster centers is then rank-ordered according 
to likelihood, and membership probabilities are used to 
mask out member galaxies in the percolation step. 

The final aspect of our algorithm is the selection of 
cluster CGs. Specifically, if the luminosity distribution 
of CGs as a function of richness is known, then one can 
evaluate the relative probability that any given galaxy 
with a color consistent with the red sequence is consistent 
with the CG distribution, with the satellite luminosity 
function; or with the foreground magnitude distribution. 



We also take into account a local density estimator, as 
CGs tend to reside in the densest parts of the cluster. In 
general, however, we find that the luminosity filter dom- 
inates in most cases. To define the filter functions, as 
with the red-sequence calibration we follow an iterative 
procedure: on our first iteration, centers are defined as 
the brightest high membership probability galaxy. These 
assigned centers are used to characterize a "first guess" 
for the luminosity distribution and local density of CGs. 
We then rerun our central selection algorithm iteratively 
until convergence. In the end, every galaxy in the clus- 
ter can be assigned a probability of being a CG, though 
typically no more than two galaxies have a significant 
probability of being the CG of the cluster. 

2.2. Cluster and Galaxy Selection 

To reduce an effect of stacking clusters of different sizes 
we choose richness cut in the catalog of 20 < A < 50, cor- 
responding to the m ass range 1.4 < M2ooc/10 14 M© < 3.6 
(|Rvkoff et al.ll2012[ ). When stacking clusters there is also 
an issue of miscentering, which we alleviate using the 
centering probability Pcg, reported in the redMaPPer 
catalog, to select clusters with Pqg > 0.95, implying we 
expect > 95% of our galaxy clusters are properly cen- 
tered. We impose no cut on redshift. 

Cluster candidate members are selected from SDSS 
DR8. To get accurate relative velocities we use only 
galaxies with spectroscopic redshift s. For each CG we do 
an initial broad selection of member galaxies by choos- 
ing all galaxies that satisfy \Av\ < v max = 4000 km s _1 
and d < 6 Mpc, where Av is the difference in velocity be- 
tween the galaxy and the CG and d is their projected sep- 
aration. These cuts are intentionally generous, since in- 
terlopers are included in our velocity distribution model. 
If a galaxy is a candidate member for several clusters, 
we assign the galaxy only to the nearest cluster, based 
on the projected distance to the CG. Finally we discard 
clusters with fewer than four spectroscopic candidates. 
To ensure that our results are insensitive to the precise 
selection criteria adopted, we repeat the data analysis 
using a large range of different criteria. We recover con- 
sistent results in all cases. 

This leaves a total of 1743 clusters containing 49,524 
candidate member galaxies (of which 13,565 have d < 2 
Mpc) that were used in the final analysis. 

2.3. Galaxy Positions Angles 

CG major axis position angles are taken from the R- 
band photometric object position angles (E of N) from 
a de Vaucouleurs model fit provided in the SDSS DR8. 
To avoid ambiguity in determining CG position angles, 
we require that the projected axis ratio of the ellipticity 
for each CG is less than 0.9. Member galaxy coordinates 
are used to determine the azimuthal angle Acf) between 
the galaxy and the projected CG major axis. Galaxies 
are put in two spatial sections according to their position 
angle relative to the major axis of the CG (see Figured]). 

2.4. Line- of- sight Galaxy-velocity Model 

Member galaxy velocities are calculated relative to the 
CG (v re \ = ^cg — ^gai)- This only makes sense if the CG 
is at rest with respect to the halo, otherwise we might 
observe a systematic offset of the mean of the velocity 
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Figure 1. Sky plot of a galaxy cluster with member galaxy can- 
didates. The CG is represented by an exaggerated red ellipse. The 
azimuthal sections used in the analysis are indicated by dashed 
lines. Member candidates are marked according to their associ- 
ated section, with green squares indicating galaxies closer to the 
CG major axis and blue triangles indicating galaxies closer to the 
minor axis. This plot contains interlopers as well as cluster mem- 
bers. 

distribution of member galaxies. If the offset is random 
the mean will be zero, but the dispersion will become sys- 
tematically larger. We find that for all stacked velocity 
distributions the mean is consistent with z ero, within ef- 
fects from gravitational redshift (see, e.g.. IWojtak et all 
([201lh . 

The velocity distribution is modeled as a Gaussian 
probability distribution with mean \i and dispersion a. 
Interlopers are assumed drawn from a flux-limited uni- 
form background, which makes the background a lin- 
ear decreasing function of recession velocity, given by 
a probability distribution with intercept n, and slope 
s(n) = (l/(2i? max ) — n)/v max fixed by the probability 
distribution normalization condition. To ensure overall 
normalization, a parameter p is introduced, which is the 
probability for a galaxy to be a member of a cluster. The 
probability for a particular relative velocity v re \ in this 
model is given by 



P(v Te i\p,a,[i,n) =p 
+ (1 



1 



■ exp 



:0rel - M) 2 /^ 



2na 

p)(s(v re \ +% ax ) +n). (1) 

Each velocity derived from SDSS includes a measure- 
ment error of the order of a few tens of km s _1 , which has 
the effect of slightly increasing the measured dispersion 
of the stacked sample of velocities. This error is small 
and should not affect the results, especially when study- 
ing azimuthal variation in dispersion. We estimate the 
likelihood for a given velocity dispersion a v by a Bayesian 
likelihood analysis, multiplying probabilities over all ve- 
locities in the data sample, 

C(p,a v ,n,n) = Y[P(vi\p,a v ,fjL,n). (2) 

i 

To generate the posterior distribution for the param- 
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Figure 2. Example velocity distribution of member galaxies along 
with the model from Equation (pQ) with highest-posterior values ob- 
tained using MultiNest (solid lines). The green crosses/line shows 
velocities along the major axis, and blue crosses/lines shows veloc- 
ities along the minor axis. 

eter a v we sample the likelihood over the prior parame- 
ter space, while marginalizing over all o ther parameters, 
using the method of nested sampling ([Feroz fc Hobsonl 
2008) . We use the publicly available software MultiNest 
(|Feroz et al.ll2009[ ). An example of a model fit to data 
for a particular radial bin is presented in Figure El 

3. RESULTS 

Figure [3] shows galaxy velocity dispersion as a func- 
tion of radius for two populations of cluster member 
galaxies; those that lie closer to the projected CG ma- 
jor axis and those that lie closer to the projected CG 
minor axis. We find a decreasing relationship between 
velocity dispersion and radius, consistent with the re- 
sult expected for a Navarro- Frenk- White (NFW) pro- 
file (jNava rro et al. 19 9^), wi th concentration parameter 
c = 5.5 (jKlvpin et al.ll2011h and typical anisotr opy of 
galaxy orbits in clusters (jWojtak fc Lokasl l2010f ) , con- 
volved with the mass function of the redMaPPer clusters 
in the given richness range. 

In Figure |4] we plot the measured difference in line 
of sight velocity dispersion as a function of radius, di- 
vided by the mean dispersion at the given radius. The 
weighted difference is calculated as Aa v / (a v ) = 2(cr™ aj — 

< in )/« aj +< in )- The dashed line (blue band) shows 
maximum likelihood (±lcr) for a constant difference in 
weighted dispersion between the two azimuthal sections. 
The posterior for the constant weighted difference is 
Aa v /(a v ) = 0.06 ± 0.02, corresponding to an absolute 
difference of A<j v = 38 km s _1 ± 13 km s _1 . This in- 
dicates that galaxies along the projected major axes of 
the CGs exhibit preferentially higher velocity dispersions 
compared to galaxies lying along the projected minor 
axes. 

We test whether the observed difference could be sys- 
tematic in nature by repeating our experiment, but ran- 
domly orienting the halos when stacking, and then split- 
ting the galaxy memberships along a random axis. As 
expected we do not detect any difference in the velocity 
dispersion of the two galaxy populations in this case. We 
also verify that our conclusions are valid if we remove the 
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Figure 3. Velocity dispersion of galaxies as a function of radius. 
Green squares are galaxies closer to the central galaxy major axis 
and blue triangles are galaxies closer to the central galaxy minor 
axis. The dashed gray line shows prediction from NFW. 

radial bin with the strongest offset (R = 0.6 Mpc), albeit 
at reduced significance of 2a. Finally, we find that vary- 
ing the cuts used to select candidate cluster members 
does not impact our results. 

3.1. Numerical Simulations 

To compare our measurement with theoretical predic- 
tions, we analyzed velocity dispersion profiles of subha- 
los in cluster-size halos from t he Bolshoi simulatio n^ of 
a standard ACDM cosmology (jKlvpin et all 1201 lh . We 
used ~ 500 dark matter halos with masses M200C > 
0.4 x 10 14 M Q and subhalos with at least 30 particles cor- 
responding to a minimum mass of M200C > 2 x 1O 9 M . 
Subhalos were identified using; Bound Density Maximum 
algorithm (jKlvpin Holtzm an 1997). For every halo we 
calculated projected velocity dispersion profiles as mea- 
sured by a distant observer. The profiles were determined 
inside two equal sections along the projected major axes 
of the halos and along the perpendicular direction, like- 
wise in the analysis of SDSS data. The major axes of 
dark matter halos were determined by diagonalization of 
the tensor of inertia computed using dark matter parti- 
cles inside the virial spheres. We considered three fixed 
orientations of the line of sight with respect to the ma- 
jor axes of the halos, with angles 5°, 20° and 90°, and a 
fourth case corresponding to a fully random angle distri- 
bution of the light cones. 

The resulting profiles are shown in Figure HI The 
dotted lines show the predicted weighted dispersion dif- 
ference, with the colored bands representing the error 
of the mean determined by bootstrapping. The mean 
anisotropy of the velocity dispersion increases strongly 
with the angle between the line of sight and the major 
axis. Moreover, it is positive (cr™ a -i > <7™ m ) at all radii. 
The signal vanishes for sight lines parallel to the ma- 
jor axis of the halo. We checked that this result does 
not depend on the alignment of the intermediate and 
minor axes of the halos. Therefore, it is a direct conse- 
quence of the fact that halos are predominantly prolate 

2 The Bolshoi simulation is publicly availab le via the MultiDark 
database (http://www.multidark.org; Riebe et al. 2011). 



(jGottlober fc Yepesll2007l ). 

Even though our calculation is based on purely dark 
matter simulation, the signal corresponding to random 
orientations of major axes of dark matter halos matches 
the observations with Xsim/dof = 10-4/8 = 1.3 (com- 
pared to xiio/dof = 15.3/8 = 1.9 for an isotropic case, 
crjf a j = a™ m ). We checked that the relative difference of 
the velocity dispersions only depends weakly on the halo 
mass and remains the same when we restrict our calcu- 
lations to the halo masses corresponding to the selected 
richness range. 

3.2. Systematics 

Miscentering of the CG has the effect of boosting the 
measured velocity dispersion of galaxies. Because this 
boost is isotropic, it will tend to dilute the signal, while 
simultaneously increasing the observed velocity disper- 
sion of the cluster stack. Furthermore, if the orientation 
of the ellipticity of the false centers is random this will 
likewise act to dilute the signal. We have attempted to 
minimize these systematics by requiring that all clusters 
be properly centered (Pcg > 0.95), but strictly speaking 
our measurement can be thought of as a lower limit. 

While CGs and host halos tend to be preferentially 
aligned in observed galaxy clusters, there is also a sig- 
nificant random component ass ociated with the angular 
separ ation of their major axes (jNiederste-Ostholt et al.l 
l2QlQh . This will tend dilute the observed signal relative 
to our predictions from numerical simulations. Thus, 
similarly in this case, one can think of our results as pro- 
viding a lower limit to the anisotropic signal relative to 
the halo's major axis. 

An additional concern is that because halo triaxiality 
depends on halo mass, the relative satellite- weighted con- 
tribution of a single halo in a stack will be different along 
the major and minor axes. In other words, the effective 
halo mass distribution for stacks along the major and 
minor axes will be different, and this could lead to an 
artificial anisotropy signal. We have estimated the level 
of this effect by use of a toy model, and found that, even 
when triaxiality scales strongly with halo mass, the effect 
is negligible. 

4. DISCUSSION AND CONCLUSION 

We present the first observational detection of 
anisotropic kinematics of galaxies in clusters of galax- 
ies. We have measured line-of-sight velocity dispersions 
of member galaxies along the major and minor axes of 
the CG, in a large sample of stacked clusters from the 
SDSS. The projected CG position angle is used to sep- 
arate galaxies into two section depending on their pro- 
jected proximity to either axes. Galaxies closer to the 
projected CG major axis are found to have a preferen- 
tially larger velocity dispersion than those that are closer 
to the minor axis, with an average difference in disper- 
sion of A<j v = 38 km s _1 ± 13 km s _1 , corresponding 
to a relative difference of Aa v /(a v ) = 6% ± 2%. These 
results are consistent with results from cosmological sim- 
ulations, with Xsim/ dof = 10.4/8 = 1.3. 

The presence of this effect is a signature of a pro- 
late velocity ellipsoid in galaxy clusters. Keeping in 
mind that veloci ty ellipsoid is a l igned with the halo 
shape ellipsoid (|Kasun fc Evrardl 120051 : iTormenl 119971 : 
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Figure 4. Difference between galaxy velocity dispersion along projected major and minor axes of the CG, divided by the mean dispersion 
at the given radius. Positive values indicate cr™ aj > cr™ m . Triangles show results from our analysis of SDSS clusters. The dashed line (blue 
band) indicates the maximum likelihood for a constant offset (zblcr errors). The dotted lines (colored bands) show results from numerical 
simulations (zblcr errors from bootstrapping), with different lines corresponding to varying selection of angle between the line of sight and 
cluster halo major axis. 



I White et al.l 120101 ) an d the halos tend to be oriented 
with the cosmic web (jAltav et al.l 120061 : iBrunino et al.l 
120071 ). the effec t is likely a remnant o f an anisotropic 
halo formation (jSparre fc Hansenll2012l : R. Wojtak et al. 
2012, in preparation). This raises an important ques- 
tion as to what extent cluster dynamics is influenced 
by the infall process and cluster location in the cos- 
mic web. The result can also be interpreted in light of 
studies indicating preferential alignment between the CG 
and its host halo (Binggeli 1982; Hashimoto et al.l 2008: 
iNiederste-Ostholt et al.ll2010[ ) , strengthening the case for 
anisotropic cluster formation. 

We furthermore find a clear correlation between the 
degree of projected anisotropy and the line-of-sight ori- 
entation of cluster halos in the simulations. This result 
opens exciting new possibilities for assessing the degree 
of random alignment of galaxy clusters. Specifically it 
might be applied to a sample of optically selected clus- 
ters to evaluate the level of selection bias. 
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